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The momentum- and frequency-dependent T-matrix operator for the scattering of electro- 
magnetic waves by a dielectric/conducting and para- or diamagnetic sphere is derived as 
a Mie-type series, and presented in a compact form emphasizing various symmetry prop- 
erties, notably the unitarity identity. This result extends to magnetic properties one previ- 
ously obtained for purely dielectric contrasts by other authors. Several situations useful to 
spatially-dispersive effective-medium approximations to one-body order are examined. Partial 
summation of the Mie series is achieved in the case of elastic scattering. 



Keywords: T-matrix; heterogeneous media; dielectric; magnetic; optical theorem; unitarity 
identity; dynamic effective medium theory; spatial dispersion. 



1. Introduction 

The transition operator, or T-matrix, of a scatterer is a basic building block of 
time-harmonic theories of single or multiple scattering It embodies information 
about the overall polarization-dependent response of a finite object, in terms of 
the incident and scattered momenta. A recent bibliographical review devoted to 
the use of T-matrices in electromagnetism (3] (mostly of the "on-shell" variety, see 
definition below) illustrates its key importance in the treatment of the response of 
heterogeneous media of various natures. Excepting scatterers of the simplest forms 
however, the T-matrix is generally a painstaking object to compute, and often 
leads to systems of equations that can be numerically problematic {e.g., Ref. 0] 
and references therein). 

In the vast majority of treatments, the T-matrix is computed on-shell. This ter- 
minology, borrowed from particle physics, indicates that the incident and scattered 
momenta have their norms fixed by the dispersion relation of the host medium. 
To be precise, the on-shell case corresponds to taking |ki| = |k2| = km with 
km = (u^/cj^EmlJ'm, where k.2 (resp. ki) stands for the arriving (resp. depart- 
ing) wavevector, u the angular frequency of the incident wave, c the velocity of 
light in vacuum, and Em and Hm are, respectively, the complex relative (frequency- 
dependent) dielectric permittivity (including conduction) and the para- or diamag- 
netic permeability in the embedding medium. These latter quantities are assumed 
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homogeneous and isotropic. With on-shell T-matrices, a scattering system can only 
be examined "from outside", i.e. by sending and receiving signals from the host 
medium However, the modern developments of multiple scattering theory {e.g., 
0, [^) have made it clear that a fully general solution of the scattering problem 
in the bulk of a heterogeneous system requires that the T-matrix be computed off 
shell, that is, with norms of the arriving and departing momenta, k2 = |k2| and 
ki = |ki| respectively, being arbitrary (not being constrained to equal km)- The 
full generality provided by the off-shell formalism allows one to study of arbitrary 
fields in bulk heterogeneous media, and also to deal with interface problems for 
which approximate techniques are available, see references in 01 • 

Moreover, off-shell computations are the proper context for discussing spatial 
non-locality (also referred to as spatial dispersion) in bulk response of heteroge- 
neous random media 0]. This allows for theoretical investigations of the various 
propagation modes that can arise as a direct consequence of the finite size of het- 
erogeneities. Information on these modes (most of them strongly attenuated) can 
be obtained either by direct computation of such "leaky modes" , or by computing 
the density of states using the imaginary part of the Green's function 

Although spatial dispersion in the electrodynamics of crystals is an old sub- 
domain of solid-state physics 0, heterogeneity-induced spatial dispersion in 
random media is a less understood matter, for which many points remain to be 
clarified |8|. More importantly, while being for a long time a pure theoretical preoc- 
cupation [9, 11-3], it has now acquired some experimental substance in acoustics 



19l |. electroniagnetism in random media [201 ] . and for electromagnetic metamateri- 
als, e.g., [H,!!!]. 

Given the complexity of off-shell T-matrices as compared to their on-shell coun- 
terparts, we restrict ourselves to the simplest three-dimensional case of a single 
sphere. Analytical results for off-shell transition o per ators in classical physics are 
the scalar T-matrix for acoustic scattering 0, 2^, 24] and the tensor electromag- 
netic T-matrix for a purely dielectric sphere 5|]. To our knowledge however, an 
explicit expression similar to that in Ref. i fo7a sphere with both dielectric and 
magnetic contrast with respect to its embedding medium has not previously been 
available in the literature. Since spatial dispersion implies the existence of an ef- 
fective magnetic-like response even in dielectric media, as has been observed by 
a number of authors [e.g., (1, H^), the dielectric and magnetic case is a neces- 
sary milestone on the road towards a realistic frequency-dependent self-consistent 
effective-medium theory. Indeedjdespite the considerable amount of work having 
addressed this issue (see Refs. [1, 0] and references therein) the latter question 
remains unsettled. 

This paper is devoted to presenting the expression of the off-shell T-matrix op- 
erator of a sphere with arbitrary dielectric and magnetic contrast, in the form of 
a Mie series expansion. This result was derived more than a decad e ag o (26| but 
remained unpublished, although having been announced in Refs. [itI. l27l] . It should 
be mentioned that by about the same time. Tip independently considered this 
same off-shell case in an abstract mathematical f;amework [H, but gave explicit 
results for a vacuum background only, and it is unclear how his result compares 
to ours. In the following, we exclude the situation of a nonzero applied constant 
magnetic field. An extension of the off-shell T-matrix formalism in this case has 
recently been put forward and exploited in connection with electromagnetic wave 
propagation in magnetochiral media (29| . 

The direct demonstration of our result can be found in Ref. [1^. However, it uses 
the approach of Refs. 0, 0] and requires a great deal of preliminary work besides 
being particularly tedious. The calculation goes in three steps. Step I: compute in 
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real space the Green's function G of the electric field in a medium containing a 
single sphere, for arbitrary positions of the source and observer points, inside or 
outside the sphere. This Green's function is obtained as a sum of four complemen- 
tary parts, each one addressing a typical situation for the emission and observation 
points, which can be independently located inside or outside the scatterer. Each 
of these parts is expanded on a basis of vector spherical harmonics (VSH) and the 
elements of this expansion involve multiplicative combinations of spherical Bessel 
or Hankel functions as is the rule with spherical scatterers. A crucial aspect of this 
calculation is that it appeals to longitudinal-electric components of the field, as in 
Ref. 0], in contrast with older approaches to the problem of wave scattering by 
spheres [s^] where these components are ignored. They have since been recognized 
as playing an important part in source regions [U], being responsible for evanes- 
cent modes originating from the scatterers. Step II consists in taking bi-variate 
Fourier transforms of this Green's function with respect to the source and observer 
positions, before applying manipulations that allow one to extract from it the off- 
shell T-matrix, see equation ([7|) where Gm denotes the dipolar Green's function in 
free space. This step is the most difficult one, since it involves non-trivial definite 
Fourier integrals on separate ranges, < r < a and a < r < oo, where a is the 
sphere radius and r is a radial coordinate. We could not express these definite inte- 
grals in closed form, but instead reduced them to a lengthy sum of explicit terms, 
added to a residual definite integral of simple form. The latter fortunately cancels 
out with an identical term arising from the VSH expansion of the Dirac singularity 
at the origin of Gm, to be subtracted from the G in the process of extracting the 
T-matrix according to equation ([7]). Hence the result can be expressed in closed 
form as in Ref. without any non-evaluated integrals. In Step III, some tedious 
re-organizations of terms are carried out to bring the result into a more usable 
form that displays all symmetry properties of interest. 

The complications makes it problematic to present a concise exposition of this 
direct approach in the case of dielectric and magnetic contrast, so it will not be pur- 
sued here. Instead, we outline hereafter a new and shorter -albeit non-deductive- 
proof of the result, which bypasses most of the difficulties of the direct approach. 
The proof consists in showing that the T-matrix satisfies a defining relation of the 
T-matrix, namely Eq. This verification only requires carrying out integrals by 
a method that can be explained relatively easily. Moreover, these integrals only 
involve simple poles whose residues can almost be read by inspection. Carrying 
out such a check remains a cumbersome task, but is a straightforward thing to do 
with a minimum amount of preliminary technical material. 

The paper is organized as follows. Our Fourier transform conventions are ex- 
plained in Appendix A. After setting up our formalism in Sec. El the coefficients of 
the Mie series of the off-shell T-matrix are given in Sec. 13.2 formulated in such a 
way that important symmetry and conservation properties 32] (among which the 
unitarity identity (33l.l33|) are made conspicuous. These properties are discussed in 
Sections 13.31 and 13.41 The principle of a proof of our result is detailed afterwards in 
Sec. 13.51 For purposes of clarity and further physical insight on the structure of the 
T-matrix, we use an intermediate decomposition of the T-matrix in intermediate 
partial "dielectric" and "magnetic" parts, which originate from similar decompo- 
sitions of the scattering potential. Before concluding in Sec. El some limiting cases 
of particular interest are examined, and new expressions relevant to applications 
to spatially-dispersive effective-medium approximations are obtained in Sec. [H 

Henceforth, the sign x stands for the three-dimensional vector product, unless 
otherwise indicated. 
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2. Green's function associated to the electric field and scattering potential 
operator 

The T-matrix is most easily expressed in the time-harmonic domain and space 
Fourier representation. The dipolar Green's function G associated with the electro- 
magnetic field in an infinite medium of relative isotropic permittivity and perme- 
ability £m and /im (the index m referring to the embedding matrix) is the retarded 
solution of the inhomogeneous wave-propagation equation [sH 



[/i-^V X V X -{uj/cfej^ G„(r|r') = \5\v- v'). 
In the Fourier representation, it reads 



(1) 



G,n(k) = ^lr. 



I -kk 



kk 



(2) 



where k"^ = {uj / cf'EmtJ'm-, the angular frequency oj being considered as a mere 
parameter hereafter. In Eq. ([1]) the permeability jjim is introduced so as to make the 
source term a pure (permeability-independent) electric current. This set-up allows 
for a consistent treatment of media with heterogeneous magnetic permeability. 

Because of translation invariance, Gm(r|r') = Gm(r — r') so that Gm(k|k') = 
Gm(k)(5(k — k'), which defines Gm(^)- In Eq. ([2]), I stands for the identity matrix, 
and Gm is expressed in terms of the transverse and longitudinal projectors with 
respect to the direction k = k/fe of the Fourier mode k. 

Consider now the one-body inhomogeneous problem in presence of a spherical 
scatterer of radius o centered at the origin. With our above convention for the 
permeability, the constitutive properties of the medium are specified by 



e(r) = em + {Es - £m)4>{'c), 



1 _ 1 /I 1 



(3) 



where Eg-, l^s are the relative permittivity and permeability of the sphere, of char- 
acteristic function 0(r) = 9{a — \r\) (9 denotes the Heaviside step function). The 
scattering potential operator U between points r and r' is 



U(r|r') = 5{r-r') 



V'x 



1 1 



(r')V' X + 



[^)\es-e^) l</.(r') 



(4) 



where the prime denotes a derivative with respect to r'; or in Fourier form [171 ]: 



U(ki|k2 



(2^)3/2 



—] ki X k2 X + (-) {Es - Em) I 



(5) 



The Green's function G associated to the electric field in the medium now obeys 
the integro-differential equation 



1 



■V X V X -{uj/cfEr 



G(r|r') = I ^^^r - r') + j d\i U(r|ri)G(ri|r'). (6) 



In formal operator notation 0], and with the help of Gm, this equation takes the 
Lippmann-Schwinger form G = Gm + GmUG. 
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3. The T-matrix 



3.1. Definitions and Derivation 



The T-matrix operator is introduced so that G = Gm+GmTGm, hence by definition 
T = U{1 — GmU)~^. This operator inversion is difficult to perform directly. Thus, 
we computed the T-matrix following Tsang and Kong 0], i.e. by first solving the 
one-body problem in real space for G(r|r') with arbitrary positions of the source 
r and of the observation point r' inside or outside the scatterer; then, by going to 
Fourier transforms; and eventually by extracting T via the relationship 

^ = Gm{G - Gm)G^. (7) 

The drawbacks of this direct approach have been recalled in the Introduction (see 
for details). Thus, a much shorter albeit non-deductive proof is provided in 
Sec. 13.51 which consists in proving that the following equation for T holds: 

T = U + UGmT. (8) 

The T-matrix is obtained as an expansion over vector spherical harmonics defined 
in Appendix A. We depart from other authors 0, Isol . 31] by using the orthonor- 



malized VSH basis {N^^, Z/„, X/„} as found in the book by Cohen- Tannoudji et 



al. |35l |. The number / > is the multipole index while — / < n < / is the angular 
number. With VL\^ as the solid angle in direction k, ^in{^\i) has the character of a 
longitudinal electric component, aligned with k, whereas Z;„ and X/„ are of trans- 
verse electric and magnetic character, respectively, and are orthogonal to k. Both 
Z;„ and are nonzero for I > 1 only. As recalled in the Introduction, longitudinal 



terms built on are necessary in any source region [3l|]. 



The following variables related to dielectric and magnetic contrast are introduced: 

Ae = {£s/em) - 1, = (nm/fJ's) - 1, = (em/es)Ae, dfi = {iis/ iim)^lJ', 

Also, we introduce suitably normalized derivatives of the Ricatti-Bessel and Ricatti- 
Hankel functions, standard in this context, defined as the product of x by the 
spherical Bessel of Hankel functions ji{x) or respectively [sH]. They read: 

- ( u \ [akajiiaK)]' (1) _ [aKhf\aK)]' , „s 

^i^a. = Vi{aK) = , . , = '{aK) = [l / 0). 

(9) 

These (p functions simplify the evaluation of our forthcoming results in the static 
limit where w — )• 0. In the above expressions, index a stands either for m, s, k, 1 or 
2 depending on the argument k^ being km = {'^ / c){emlJ'mY^'^ , kg = {uj/c){£sfJ.sy^'^ , 
k, or the outgoing or incoming momenta ki or k2 respectively. Finally, let 

c, _V^i,a-'Pi,p p _ - k^^^i,a _ ii{aki) ji{ak2) 
^''"^= kl-kl ' ^''"^ = kl-k} ' '^''^^ = ^fc^ 

Useful limiting behaviors are ^i{x) = (/ + 1) — x^/ {21 + 3) + O(x^), so that Si^ap ~ 

-c?l{2l + 3) and Ri^fi \ when a ^ 0, and ^f\x) = -I + x^/{2l - 1) + 

0(x4) + 0((ia;)2'+i)|l] 



^ For this reason, slightly different normalizations, in the form of alternative functions Qi{x) = ipi{x) / (l + l) 
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3.2. Off-shell components 

Our VSH expansion of the T-matrix, which is the main result of this paper (see 
Introduction), reads 



T(ki|k2) = T/'^^ihlk^) A,„(f]kjBr„(J7kJ, (11) 



i>0 
A,B = N.Z,X 



n=~l 



where Ti^^(ki\k2) = T{^^ {ki\k2) = ^(A;i|A;2) = Ti^^{ki\k2) = due to spherical 
symmetry, and 



3 ;„2 



2a-^ k 



T[''\ki\k2) = —^5£ 



l{l + l)5e 



_{^m/£s)'^l,s - 



(1) 



+ Ri.i2 - 1 



Tf'^{ki\k2) ^ rf^(fc2|A;i) ^ 2a3^ 
^/liTTT) ~ + 1) ~ ^ /"m 



5e 



{£m/£s)(Pl,s - f i 



(1) 



+ 1 



Jl 



(12a) 
(12b) 



Tl^^{kl\k2) 



2a^ fcm 

vr firn. 



(SeRi^u + 5iJLklSi^is){5£Ri^2s + SfJ'k2Si,2s) 
kfk2 {kf — k'^)Si^is — (k'2 ~ k'^^)Si^2s 



+ 6n 



U2 



i^l '^2 



+ 



5e k^ikf - ki)Ri^u - kfiki - ki)Ri,2s 



U2 

A/™ 



1^1 "-2 



Ji,i2, (12c) 



Tl''''{k,\k2) 



2ar kik2 

71" A'm 



{l^llRl^ls + A:^A£5;,l.)(A^i?;,2s + fc^A£5;,2s) 
{Hrn/fJ^s)^l,s - 

(A/ii?/,i2 + k^AeSi 



(12d) 



(AMfc? - e A£)(AMi - fc^A£)%-^^^ - - ' 



i-2 
"■l "'2 



12- 



Elements T^^"^ are defined for I > 0, whereas the other types are defined for I > 1. 
At the price of additional algebraic manipulations, we checked that in absence 
of magnetic contrast (fig = fJ-m), these compact and symmetric expressions are 
equivalent to those by Tsang and Kong 0] (who use non-orthonormalized VSH). 
The equality of Tj^^{ki\k2) and T^^^ {k2\ki) is a consequence of reciprocity [3l|. 



3.3. Basic properties 

Briefly, the main properties enjoyed by these matrix elements are as follows. First, 
j'NN ^ rpzz T(^-^ are symmetric under interchange of ki and k2, and are non- 
singular for all finite values of {ki,k2) (possibly complex). With respect to this 
property, note that at a zero of ji{aka), a = 1,2, the product ^pi^ajiio-ka) is always 
finite. In addition, the triple limit /ci — )• 0, A;2 — )• 0, a; — >• is uniquely defined, 
because the limits commute. Also, the limit 

Urn hm T^^^ {k,\k2) = lim T^^{k\k), (13) 

Kl— >Km K2^Km k^k,„ 



and Q^j^\x) = —(p^^\x)/l were used by us in Ref. [17l |. 
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uniquely defines the so-called 'on-shell' elements (a concept relevant to transverse 
components only, see below). 

The denominators identify the NN, NZ, ZN and ZZ terms as electric multipole 
contributions, and the XX terms as magnetic multipoles. The prefactor kik2 in 
T^^^ {ki\k2) is results from the operator ki x in the potential expressed by 
Eq. ([5]). As in classical Mie scattering (i.e., for the on-shell T-matrix, see Sec. 14. 2p . 
transverse electric and magnetic polariton resonances l3, 13, 3^ arise at complex 



frequencies for which denominators vanish, namely when: 



^™ ^liaks) - \ ipl'\akm) = {l> 1). (14) 



3.4. Unitarity identity as a consistency check 

In the non-dissipative case, all constitutive parameters and U(r|r') are real. Imag- 
inary parts of the T-matrix elements, rooted in the outgoing- wave prescription 
+^0+ in Eq. ^ that defines Gm(k), arise solely from the V'z^m- The fact that these 
imaginary parts are separable in the momenta ki and k2, see Eq s. (|12ap - ()12dp . 



is deeply connected with the well-known unitarity identity j32l434i| . This identity 
is a generalization of a well-known statement of energy conservation in scattering 
theory, [i^ , wherein the scattering as cross-section is equal to the extinction cross- 
section of the scatterer cje in the absence of absorption. As relationships (jl7p below 
show, this identity is deeply connected with the fact that the sum of the second 
and third terms enclosed in braces in Eq. (fHcP or (fTMl) reduce to the left or right 
factor in the numerator of the first term in the same expressions, when ki = k^ 
or k2 = km, respectively. This provides an easy consistency check for expressions 
in Eqs. (jl2ap - ()12dp . For definiteness and further reference, this unitarity identity 
is derived in Appendix ICl in the vector case. Its generic operator form reads 



l(r-Tt) =rtlm(G^)r = TIm(G^)rt. (15) 



For spherical scatterers, it takes the form (see Appendix ICl) 



ImJ {ki\k2) = '^Hmkm y"dOqT(ki|A;mq)(l -qq)T*(/cmq|k2). (16) 



Ensuing identities for matrix components are obtained as follows. Since for a sphere 
= T^'^ = = T-^^ = 0, the sums over n in Eq. (jlip are purely real 

(see AppendixiBjfor their explicit value). Then T*(ki|k2) is expressed by Eq. (fTTj) 
provided that Tj^^ {ki\k2) is replaced by Tj^^* {ki\k2) in this expression. Then, ex- 
panding identity in Eq. (I16p on the VSH basis and identifying mutually orthogonal 
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components leads to the following relations, to be obeyed for each /: 



Imr;^-^(A;i 


k2) 


~ '^f^mkm-L I [kl 


\km)J-l {k'm\k2), 


(17a) 


'7 '7 , 


k2) 


— '^i^mk-mXl \k\ 


Z Z^ / 7 

km)Ti {km\k2)-, 


(17b) 


Imrf^(fci 


k2) 


~ '^l^mkmXl {kl 


kM^''*{k^\k2), 


(17c) 


Imr,^^(fci 


k2) 


t, k T^^ik. 


km)T^^ {km\k2)i 


(17d) 


Imr,^^(fci 


k2) 


= '^fJ'mkmTi {kl 


km)Ti {km]k2)- 


(17e) 



For I = 0, the last equation must be replaced by ImTg (A;i|A;2) = 0. These relations 
can be explicitly checked on the matrix elements themselves with the help of the 
formulas (the second equality stems from the Wronskian [s^ W{ji,hi^^) = i/x"^): 

lmip\^\x) = [x\h\^\x)\'^]~\ (18a) 
^i{x) - v?p)(x) = [ixh\^\x)ji{x)]~\ (18b) 

We close this section with the following remark. Combined with the symme- 
try properties of the matrix elements, Eqs. ()17ap . (Il7b|) and (I17ep imply, for 
real k and non-dissipative media, the positivity of ImT;^^(A;|A;), IvaT^^ {k\k) and 
ImT;'''""'^(fc|A;). Moreover, setting t;(/ci|/c2) = (7r/2)/im^m^K^i 1^2) where Ti stands 
for either T^^ or T^-^ , it is easily seen that identities in Eqs. ()17aP and (|17bp 
imply the existence of real symmetric functions ti{ki\k2) such that 

n{ki\k2) = U{ki\k2) + 2iJlihjMlMt^, \ti{k^\k^)\ < (19) 

1 + ^l-At^{k„^\kra) ^ 

a choice of sign in front of the square root having been made. In words, the real 
part of T^^ or T^^ fully det ermines the latter quantities as functions of ki and 
k2- The w-dependence of T is closely tied to its fc-dependence, as shown by the 
way km occurs in the imaginary part of Eq. (jl9p . Finally, the inequality in Eq. 
(I19p allows one to define phase shifts bi (real in absence of dissipation) different for 
the ZZ and XX components, such that Ti(kra\km) = sin((5;) exp(i(5;) [j^. Isil. [4o| . 
Both ti and bi are analytically continued as functions of e^j ^m, P^s and [im in the 
dissipative case. 



3.5. Sketch of a proof of Eqs. (T^ 



We prove Eqs. (I12p by verifying relationship ([8]), taken as a definition of T. To pro- 
ceed, and anticipating further applications to heterogeneous media with spherical 
inclusions, it is convenient to split up the scattering potential U into its 'dielectric' 
and 'magnetic' parts, f/^ and U^, defined from Eq. ([5]) by alternatively suppressing 



the magnetic or the dielectric contrast: = U\^^=fj,^ and = U\s^=err, |27|]. Our 
aim is to show explicitly that 

T={U' + U'GmT) + ([/'' + U^'GmT), (20) 



which is equivalent to proving Eq. ([8]). 
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From this perspective, VSH representations of C/^ and are obtained from 
the T-matrix components of Eqs. (jl2p by keeping only their lowest-order term 
in an expansion in powers of the dielectric and magnetic contrasts Ae and Afj,^ 
since U^''^ are proportional to these quantities. The following nonzero elements are 
obtained: 



U[''^{h\k2) = U[^''{k2\h) 



'^yI(^Tl)^AeJ,,l2, 



TT 



^e-Ki,12 Jl,l2, 

TT IXrn 



Ul'^\k^\k2) 



U['''''ik,\k2) 



2a-^ kik2 2 A c T 

A^5/^12 Jz,12, 

2a^ kik2 

IT flm 



AfiRl 12 Jl12- 



(21a) 
(21b) 
(21c) 
(21d) 
(21e) 
(21f) 



Next, one must compute U^GmT and U^GmT. The calculation goes as follows: 
first, expand for instance / U^(ki|qi)Gm(qi|q2)T(q2|k2) d^gi d^g2 on the VSH basis. 
Using the orthonormalization properties, we arrive at expressions such as 



{U'GmT)f^{ki\k2) = 4vr / dqq' g' {q)Uf ^ ^ {k^\q)Ti^ ^ {q\k2) 



teZZ, 



+ g''{q)Uf^''{kMTr{q\k2 



■NZ, 



where g and g are the transverse and longitudinal parts of that we have 
written Gm(k) = g'^ {k){ \ — kk) -|- 5f^(A;)kk. Such expressions only involve prod- 
ucts of the form g^ {q)U[ {ki\q)Ti^^ {q\k2), or g^ iq)U[ {ki\q)T^^ ^ {q\k2) or 
g^{q)U[^^{ki\q)Ti^^{q\k2), where A,B = N,Z,X. Upon going back to defini- 
tion Q of the function (pi{x), one observes that these integrals over q all reduce to 
generic contributions of the type 



I = ^"dg/(g2)|a(g2)j,(ag) + /3(g2)[agi,(ag)]'}{7(g2)j,(ag) + 5(g2)[«gj-,(ag)]/|^ 

where /, a, /3, 7 and 6 are rational functions of q^. These integrals are compute d by 
the following standard method in presence of trigonometric or Bessel functions [4ll | . 
Splitting up I into two equal parts by writing / = 1/2 + 1/2, then performing al- 
ternatively the substitution ji{aq) = [/i^^ (aq) + /ip'' {aq)]/2 in the right-hand factor 
(in the first instance 1/2) and in the left-hand one (in the second instance), and 
appealing next to the change of variable q — )■ —q, with hf'\—aq) = (— l)'/ip^(ag). 



^This requires Taylor-expanding 5e and <5/i in powers of these quantities, too. 
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yields the equivalent form 



I=\ Mf{q%{a3l{aq) + n<^q3l{aq)\']{lhf\aq) + 5[aqhf\aq)]'] 

+ [ahf\aq) + P[aqhf'^ (aq) ]'} {7 ji{aq) + 6 [aq jiiaq)]']). (22) 



To be precise, we indicate that this transformation turns products Ri^iqji{aq) and 
Si,iqji{aq) in the original integral into the following new quantities: 



2 r 2 (1) 

RiM^Mq) 2 U — "-i (23a) 

q — K-^ 

(1) _ 

SimM^Q) ^ ^^-^^S'V'?)- (23b) 
q — 



Once cast in the form of Eq. (|22p the integral can be computed by contour inte- 
gration, closing the integration path on the real axis using a half-circle of infinite 
radius in the upper half-plane. Since the transformation generates products of ji 
and functions, the contribution of this half-circle vanishes. Besides poles 
due to the transverse part of the Green's function (if present), the transformation 

endows the integrand with a single pole at g = due to products ji{aq)h!p {aq) , 
and double poles among ifes, zt/ci or ib/c2 because the functions R and S have 
been modified according to Eqs. (|23ap . (I23bp . The pole (7 = must be handled 
by a principal value prescription, whereas in the pairs of poles of opposite sign 
that of minus (resp., plus) sign is shifted in the lower (resp., upper) half-plane by 
an infinitesimal amount. In this way only poles with plus sign contribute. Most 
often in these integrals the associated residues can be read by inspection. Exten- 
sive use is made of the Wronskian identity ()18bp in subsequent reorganizations to 
reduce residue contributions coming from terms such as Eqs. (j23ap or ()23bp . For 
instance, the pole q = ki generated by ()23bp gives rise to a residue proportional to 

{(fl^i — ip[^^)h^^\aki), equal to l/[iakiji{aki)] by virtue of Eq. (Il8bp . In general, 
the remaining function ji{aki) in this denominator cancels with a similar factor 
present in the numerator of the multiplying term within braces in the integrand of 
Eq. (|22p . in which the ji have not been transformed. Albeit lengthy, the calculation 
is thus straightforward. 
Adding, respectively, the contributions of potentials U^, read from Eq. (|2ip 
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to the matrix elements of U^GmT and U^GmT computed by this procedure yields: 



(24a) 



{U' + WG^T)f\ki\k2) = {U' + WG^T)t'\k2\k{) = Tl'^{ki\k2), (24b) 



{W + U'G^T)f\ki\k2 



3 1.2 



2a^ fe, 



5eRi^is{5eRi^2s + 5nk'^Si^2s) 



(1) 



+ 



Se k2(ki kj^)Ri^is ^i(^2 ^m)-^l, 



2s 



U2 



J/,12, (24c) 



{U' + U'G^T)f''{k^\k 



2a? kik2 

vr Urn 



kl^AeSi,u{AnRi,2s + kl,AeSi,2s) 



(1) 



klAe{A^^kl - fc^A.)%-^ - k^AeS^u 



9„3 1.2 

([/M + u^'G„,T)f^{ki\k2) = —2211 



J; 



/,12i 



(24d) 



6fikfSi^is{6eR^2s + '5/^^25'z,2s) 



(1) 



^1^2 (^1 f^m)^l,ls (^2 ^m)^l,2s 
T-2 



(C/- + [7-G,„T)f^(fci|fc2) = — ^ 



A^iRi^u{A^Ri,2s + kl^AeSi^2s 



Ji,i2, (24e) 



(1) 



Aiik^{Ai2k2 - fc^Ae)— 2 -2— 



^U2, (24f) 



other matrix elements being zero. Identity (j20p can now be checked by mere in- 
spection by comparing these expressions to the matrix elements in Eqs. (jl2p . 



4. Limits and values of interest 



Some particular limits and values of interest are now examined. The sphere volume 
is u = (47r/3)a^ 



4.1. Point-like limit 

In the mathematical "point-like" limit where the sphere radius a goes to zero, the 
matrix elements in (1121) reduce to 



(2vr)3 fJ-s + 2/i, 



W (^il^2) = 7TTZV^^ :\ r ^n- (25b) 
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and T, 



2T, 



I 

•NN pt. 



NZpt 



{ki\k2 



V2T[^^^\ki\k2) while rf^P*(A;i|yfc2 



; {ki\k2). To lowest order in the sphere radius, the T-matrix thus reads 

1 



TP*(ki|k2) 



(27r)3 



{uj/cf£mas\ - — ki X k2X 



(26) 



where a^, and are the quasi-static electric and magnetic polarizabilities of a 
sphere [361]: 



Of: 



47ra 



3 



/is + 2/i, 



(27) 



m 



The quasi-static expression for electric polarizability is too crude to obey the uni- 
tarity identity, and a variety of prescriptions have been developed in recent years 
to include finite frequency corrections to the point-like model that satisfy unitarity 



42l |. Corrections to the quasi-static limit that satisfy both unitarity and causality 



were developed in Ref. [17 1. 



4.2. Transverse on-shell elements 

For scatterers immersed in a homogenous background media, the calculation of 
physical quantities proceeds via T-matrices sandwiched between the homogeneous 
media Green's function Gm (s-g-, GmTGm), the poles of which select the "on-shell" 
T-matrix elements with ki = k2 = km- In the scattering and extinction cross- 
section calculations of Eq. ()C6|) . the on-shell T-matrix elements are proportional 
to the Mie coefficients classically obtained by solving the exterior problem where 
the source and the observer both lie outside the sphere [s^]- Specialization to 
this case of expressions (jl2cp . ()12dp after a few reorganizations that involve the 
Wronskian identity ()18b|) . yields the standard values of these coefficients, which in 
our notations reads (see also [l7| ) 

i^'u k T^^ik \k \ - -^'^"^"^^ ^m.m-eram.s , . 

l'^ H-mi^m^ I \r<"m\'^m) — — TJ\ {^OAJ 

^ h} '{akm) fJ.m'^Ls - fJ-s^Pi^m 

The right hand sides of these equations are the dimensionless T-matrix elements 
typically manipulated in on-shell theories. The factor 7r/Xm^m/2 arises from slightly 
different conventions and normalizations that are generally practiced between off 
and on-shell theories. 



4.3. Equal momenta 

The case of forward scattering ki = k2 = k is particularly important for applica- 
tions to random media, since it is the one relevant to the computation of the first 
correction in the volume density of scatterers, to the non-local effective permittiv- 
ity and permeability of the medium [1, [1^. In this case, the Mie series that defines 
T can be partially re-summed. Only the result is presented here, the calculation is 
carried out in Appendix [Dl 
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Symmetry considerations allow us to decompose T(k|k) into longitudinal and 
transverse parts as 



T(k|k) =r^(A;)kk + T^(A:)(l-kk), 



(29) 



where T'^{k) = {k) + T^{k), and where T^, and T-^ are defined by: 



5^T,^^(A^i|fc2)5^N;„(fikjNr„(OkJ = r^(A.)kk, (30a) 

l>0 n 

Y,Ti^^{k\k)Y,^inm^Lm = r^(A;)(l -kk), (30b) 

1>1 n 

5^T;^^(A;|fc) j;X,„(J]k)X;„(Ok) = r^(A;)(l - kk). (30c) 



The sums over n are given by formulas (IBSp . from which we deduce that: 



^'(^) = E^^^'''''(^i^)' T{'ik) = Y.^^^Ty^^'m (31 

l>0 1>1 



where the matrix elements, read in Eqs. (jl2p at unequal momenta, can by evaluated 
in the limit /ci , A;2 — )• /c by means of Eqs. (|D5p . Introduce now the function 



5(x) = -^(2/ + l)(^Ka 



, l-Jo(2x) 
2x2 



(32) 



which is such that S{x) = 1 — ^rr^ + 0{x'^). Appendix B shows how part of the 
sums over / that result from the above limiting process can be expressed using 5. 
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One ends up with: 



V 



(2 



L>i "pi^s - {es/smmX 



1 2 



ak 



+ 1 



(33a) 



^2 



6e + 



s {Se-6fx)^k^ 



S{ak) 



+ 



H'Tn 



c- (1.2 l.2\2 



jijak) 
ak 



(33b) 

2 



(2vr) 



-r^(A:) = — 



+ ^(^/^^' - ^^f^s) (^^) - 1] (33c) 



+ 



(A//A;2 - fe^Ae)^ 
(A:2 - A:2)2 



ak 



Expression (|33ap is the frequency-dependent counterpart of the static momentum- 
dependent expression obtained by Diener and Kaseberg [ij], to which it reduces 
when cj — )• (see also Eq. (48) of Ref . |^ ) . Apart from the occurrence of different 
magnetic permeabihties in km and ks that enter the definitions of 99 and ipi^g, 
this longitudinal term has the same form as in the case with no magnetic contrast. 
In Ref. |8|, the expression of T^(/e) provided in the case = /^s involves integrals 
that are left unevaluated. Instead, the present result is fully explicit: the transverse 
part at fis = Mm follows from using this equality and setting 5/i = A/i = 0in the 
above expressions. 

Though this is not obvious from the above, expressions of {k) and T^{k) are 
regular in the limit k ^ kg. This can be shown by using Taylor expansions, more 
particularly expansion (jPSp . In this case, it is actually easier to check regularity 
term-by-term in each individual term of the non-resummed Mie series, see Eqs. 
(|D6al) and (lD6bj) . to which one can always go back in case of problems in numerical 
evaluations near this limit. 

It should finally be noted that the right-hand side of Eq. (I33ap goes to infinity 
in the limit where — )• (which is almost the case at the plasma frequency in the 
high-frequency limit of dielectric response [1^), unless k = Q. Then indeed 



T'^ik) ~ 



V k^ Eyn 



(27r)3 ^m £s 
3 



3^/(2/ + 1) 



i>i 



jijak) 
ak 



[2ak Si{2ak) + jo{2ak) + cos{2ak) - 2] - 1 L (34) 



(27r)3 nm Es l4(a/c)2 
and the function within braces, which arises from formulas taken from (l^ and 
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where Si(x) is the sine-integral function, has no other real zero than A; = near 
which it behaves as — (afc)^/15. 

5. Concluding remarks 

We derived the off-shell T-matrix of a dielectric and magnetic sphere, provided 
relatively simple means to check this result, and some particular limits of physical 
importance were examined. Leaving applications to further work, we close with the 
following remarks. 

First, the introduction of the intermediate functions <pi, Ri and Si was found to 
be a quite useful device in trying to put some order in the structure of our results, 
and in helping displaying physical symmetries of interest. 

Second, it is observed that only the magnetic extension allows one to recover 
the limiting case of perfectly conducting inclusions: as is explained in Ref. [13] (p. 
790), this ideal case corresponds to formally taking the joint limit — )• oo and 
/Xs ^ in the scatterer. Such limiting values allow one to retrieve from Eqs. (lC6b[) . 
and (|28ap . ()28bp the well-known Mie-Debye low- frequency scattering cross-section 
of a perfectly conducting sphere obtained from Leontovich's boundary condition 
with surface impedance Z = {e.g., Ref. [sH], formula 16.159). This cross-section 
is larger by a factor 1.25 than that found for £^ = 00 but fig = /x^- Similar limits 
can easily be taken in the off-shell expressions. 

It should be remarked that even though the most useful physical quantities are 
obtained from either on-shell matrix elements, ki = k2 = km (e.g., cross sections) 
or forward scattering, ki = k2 = k, for effective-medium approaches in random 
media, it was only by computing first the T-matrix at unequal momenta that 
we can currently reach in explicit form these quantities of interest. This should 
be clear from the definition T = U + UGmT, that involves an integration over 
arbitrary momenta. In the purely dielectric case, an alternative method has been 
recently proposed @] to directly obtain the relevant elements at equal momenta, 
but the outcome involves integrals to be done numerically, and the method has not 
yet been extended to magnetic contrast. In this respect, an appealing perspective 
might consist in comparing our results in absence of magnetic contrast to that 
of Ref. @] to the purpose of deriving identities for these integrals. This might 
ultimately lead to a shorter path to obtaining T-matrices at equal momenta in 
other cases of interest beyond dia- or paramagnetism. 

Finally, the behavior in the limit — ^ emphasized at the end of the previous 
section indicates that in this case for k ^ 0, the perturbative approach that consists 
in computing the effective longitudinal dispersion relation of a composite medium 
to one-body order 0] would fail, since the longitudinal part of the T-matrix goes 
to infinity. Singularities also arise at polariton resonances. In situations of the sort, 
it has sometimes been found that in effective constitutive parameters, the first 
correction to the homogeneous matrix changes its usual proportionality to /, the 
volume fraction of inclusions, into a proportionality to some lesser power of / (e.g., 
[3]). Such cases therefore deserve special attention when considering applications 
of the present results to effective-medium theories. 

Appendix A. Fourier transform conventions 

Our Fourier transform conventions are as follows. This work makes use of generic 
operators, say A(r|r'), which may contain derivatives, with "input" point r' and 
output point r. By convention, their space Fourier transform is taken up by mul- 
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tiplying on the right by a factor e~*'' ''/(27r)^/^, and on the left by e^"^^''^' , 
and by carrying out the integrals over r and r' in the infinite volume to obtain 
the transform yl(k|k'). This use of a normalized plane- wave basis is standard when 
dealing with operators. 

However, whenever A(r|r') = A[y — v') is translation invariant (we use the same 
A by abuse of notation), we write A(k|k') = (5(k — k')^(k), which follows from 
computing A(k) as the transform of the one-entry function ^(r) by multiplying 
the latter by e~^^'^ and by integrating over r. This is the standard practice of 
solid-state physics. 

In the present context, this use of two conventions is necessary to spare us from 
dragging factors (27r)^/^ in translation-invariant expressions of interest expressed as 
Fourier transforms. No confusion will result since the use of the operator convention 
is indicated by the vertical bar between two variables. 



Appendix B. Vector spherical harmonics 

The Vector Spherical Harmonics used in this work are defined for / > and — / < 
n < I as [3f 



N;„(Ok) = ky;„(f^k), (Bla 
1 

1 



Z;„(Ok) = -j=^==VnM^^). (Bib) 
X,„(f]k) = rr^ X VnM^^)- (Blc) 



where the Yj„(0) are the usual scalar spherical harmonics j36|], and where V^^ is 
the angular part of the differential operator V = \s.d/dk + (l/A;)Vf7^ in spherical 
coordinates. Another standard notation for the VSHs is Y^^^''-^^'^'"^ (see e.g., [3^). 
However, the present notation, already employed by us in Ref. [13], alleviates the 
need for superscripts. 

The VSH are such that Xi„(r2k) = kxZi„(Ok) and Z,„(J7k) = -kxX,.„(17k)- Ob- 
serve that Xqo and Zqo are identically zero. Under parity, X;„(— k) = (— l)'X;„(k), 
Ni„(-k) = (-l)'-iN,„(k) and Zi„(-k) = (-l)'-iZz„(k). These VHS are or- 
thonormalized: 

y"dQkAL(J7k).B;.n'(f^k) = 5a ,3^1,1' 5n,n' ■, (B2) 
where A, B stand indifferently for N, X or Z. The closure relationship reads: 
^{N,„(Oi)Nr„(02) + Z,„(l^i)Zr„(1^2) + X,„(l^i)Xr,(1^2)} = \5{^i-VL2). (B3) 

In 



Introducing u = ki-k2 and the Legendre polynomial Pi{x) defined by the generating 
function (1 — 2tx + = X];>o Pl{^)^\ the following sums are obtained {e.g., 
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0): 



2/ + 1 

J^N;„(OkjNrjl]kJ = ^_P,(u)kik2, (B4a) 

n 

5;N;„(J7kjZrjl7kJ = ^lli-^ P/(n)ki(ki-nk2), (B4b) 
V 47rV/(/ + 1) 

2Z + 1 

Z,„(J]kjNrj(7kJ = — ^===P/(n)(k2 - ^ki)k2, (B4c) 
^ 47rV/(/ + 1) 

2/ + 1 

X] Z;„(J7kjZr„(f)kJ = 4^;(;^i) -PrWfe - ^xkl)(kl - ^£2) 

2/ -1-1 -V- 
+ 4^;(; + i) -P/(^)(l - kiki - kska + ^/kika), (B4d) 

X^X;„(QkjX?JQkJ = X k2)(k2 X ki) 

+ 4|m)^'^^)^^'-^^^^)- 

Since = 1 and -P/(l) = Z(Z + l)/2, the only non-zero sums at equal angles 

S^ki = f^ka are: 

5^N^„NL = ^kk, (B5a) 

n 

2/ _|_ 1 

j;Z,„Zr„ = ^Xz„XL = ^(1 - kk). (B5b) 



Appendix C. Unitarity identity 

Let A' the complex conjugate of operator A; A^ its transpose; and A'^ its Hermitian 
conjugate, in the direct or Fourier representations: Since A-{yi\y2) = ^jj(ri|r2), 

^^^.(ki|k2) = ^ji(ki|k2), A\.{vi\v2) = A*..(r2|ri), and 4j(ki|k2) = ^*,(k2|ki), 
operators and ^ commute with Fourier transforms. Our first step is to express 
by a condition on U the reality of the constitutive parameters. Potential ?7, as 
a generalized response function, is subject to Onsager's symmetry principle for 
kinetic coefficients that translates here into the principle of inverse propagation 
of light (or reciprocity). Assuming the absence of a constant external magnetic 
field, this reads: U(ri|r2) = U (r2|ri), or U(ki|k2) = U (— k2| — ki). Meanwhile, 
absence of dissipation translates as U(ri|r2) = U'(ri|r2), or U(ki|k2) = U'(— ki| — 
k2). Combining both sets of equalities implies that U = U'^ . Therefore, G obeys 
- = Uy = ul = G\^^ - Gt-i, so that using G = G^ + GmTGm provides: 

G^{T + TGlT^)Gl = Gm{T^ + TGmT^)Gl. (CI) 



The desired unitarity identity on the T-matrix follows 33|, [34 1 : T — T^ = {Gm — 
gI)T = T{Gm-Gl)T^. Noticing that gI = G*^, one ends up with equation (fT5|) . 
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With ([2|), we have for real km'- 

ImG™(k) = |^(|-kk)5(A;- fe^), 

ImG™(r) = ^ / dOk( I - kk )e^'-^-\ (C2a) 

Since Ty also obeys Onsager's principle, relations (fTSl) can be rewritten as 

ImT(ri|r2) = y"d3xi As T(ri|xi)Im (Gm(xi - xs)) T*(x2|r2), (C3a) 
l[T(ki|k2)-r(-ki|-k2)] (C3b) 
= lj^^J'mkm j dr2gT(ki|A;„q)( I - qq)T*(-A;„q| - ks). 

Scatterers for which the origin of coordinates is a symmetry center obey the prop- 
erty, inherited from [/, that T(— ki| — ks) = T(ki|k2). Equation ()C3P then entails 
Eq. (jl6p in the main text. 

To retrieve the unitarity relations in their usual form, let the incident field be of 
the form Ei(r) = e*^'''''"~'^*^Ei, with |ki| = km- The scattered field at large distances 
from the scatterer, Eg, is such that: 

E frl /■ (\\d'^V i{k„,r-uit) 

(C4) 

where kj = /cm?- The total field reads Etot(r) = Ei(r) + Es(r). Accordingly, the 
total complex Poynting vector is the sum of its incident, scattering, and extinction 
parts: Stot(r) = Etot(r) xH*„t(r) = Si(r)+Ss(r)+Se(r), with Si(r) = Ei(r)xH*(r), 
Ss(r) = Es(r) X H*(r), and Se(r) = Es(r) x H*(r) + Ei(r) x H*(r). Denoting 
the time-average of the real incident Poynting vector by (Si)(r) = ^ReSi(r), the 
scattering and extinction cross-sections are respectively 

where surface integrals a re p erformed on a sphere whose radius goes to infinity. 



centered on the scatterer [36|]. The outer medium being lossless, we find after some 
algebra that 

tTs = ^^E*- /df^/Tt(ki|k/)(l-k;k;)T(k;|ki)Ei, (C6a) 
= i^!^§JliE* • i [T(k,|k,) - Tt(k,|k,)l %- (C6b) 



kmEi 



Hence Eq. (jlSp implies the weaker conservation statement = as where the T- 
matrix is evaluated on-shell with ki = kf = km- 
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Part of the terms in the Mie series of the T-matrix can be exphcitly summed. 
These are terms with no exphcit frequency dependence. The calculation consists in 
identifying and reducing them, appealing to well-known sums that involve spherical 
Bessel functions to produce closed-form expressions. 



D.l. Longitudinal part 

Observe that the longitudinal part involves the frequency-independent terms 



E 

l>0 



1^1 "-2 



I aki ak2 



n=—l 
I 



^ kiji{aki)ji'{ak2) - k2ji{ak2)ji'iaki) ^ w* rn \C C m^\ 



l>0 



Since 



31 {aki )ji '{ak2)Yir, (Ok, )Yil{^^. 

In 

1 d 



1 d 
(47r)2 d{ak2) 

k2 - A:iki.k2 



A ^( I Nio(a|ki - kal) = -— ji(a|ki - kal) 
A-K o{ak2) 47r |ki — k2 



■(ki.k2)kik2. 



(D2) 



the sum in Eq. (jPlh evaluates to 

1 ji(a|ki - k2| 



47r a|ki — k2| 



(D3) 



With u = ki.k2 and on account of Eq. ()B4ap . this longitudinal part reduces to 
T^^(ki|k2)^ j;r,^^(A;i|fc2) Nz„(J7kjNrjf)kJ (D4) 



Z>0 



n=—l...l 



(27r)3 



l{2l + l)5ePi{u) ji{aki) ji{ak2) , ji(a|ki - k2 



+ 



alki - k2 



kik 



l«-2, 



whence expression in Eq. (|33ap for ki = k2 = k. 



D.2. Transverse part 

For arbitrary ki and k2, contributions involving VSHs Z and X do not simplify as 
easily. Still, for ki = k2 = k, some partial evaluations of contributions to the Mie 
series are possible. Let us first write down a few useful limits, making use of the 
following derivatives: 

dSi^ks _ ^ f c dRi^ks _ kl / dipi^k \ 

dk^ ~ k^- kl V dk^ ~ k^- kl V dk^ J ■ 
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lim 6/ 12 = ^To-, 



lim i?;.i2 



lim 



[kf — k^)Siis (^^2 k'^)Si2s 



lim 

ki ,k2^k 



'"'1 '^'2 

S,,i2 + {kl - kl) 



Lis 



L2s 



A-2 _ 1.2 



k'^-kl 



Slks 



"1 



(D5a) 
(D5b) 



^^l^k .,2 _ ,2 

(D5c) 



lim 

ki,k2—^k 



^2(^1 ^m)^',ls ^1(^2 k'^^)Ri2s 



^2 t-2 



lim 



At]^ "-2 



2,4 ^ ^l,k , /,2 ,„2 x,.4 d Rl,ks 



' dk^ B 



+ {kg k^)k 



k-'Rik, + ek^ 



Si 



ks 



dk'^ k"^ 

d(pi,k 
dk"^ 



(D5d) 



The above limits allow us to write the transverse elements at equal momenta as: 



Ti^^{k\k) = 



la-^ki {deRi^ks + Snk'Si^ksY 



'A r 



+ 



1.2 



(.kt k^)—^ (be 



dk^ 



ks 



dip, 



Lk 



{fJ.m/fJ-s)^l,s - f i 



(1) 



{^^,k^ - klAe?^-^ + ^f^k' ^ ''''' 



9A;2 A;2 



5fc2 



(D6a) 



1 


' ji{ak)' 




ak 



(D6b) 



ii{ak) 
ak 



The second and third terms of both these expressions depend on / only via ipi^k^ 
d^i,k/dk'^ and ipi^s- The last step consists in appealing to S{x) defined in Eq. (j32|) . 
and to the following result: 

Y,{21 + l)^[Mx)? = I [Six) - 1] = -Ix^ + O (x^) , (D7) 
i>i ^ 



whereas the evaluation of X^i>i (2Z+1)(/Jz,s [ji (afe)]^ by means of elementary functions 
is most probably not feasible (this function should admit, for all I, all the zeros of 
ji{aks) as poles relatively to the variable ks). Using S{x) and Eq. (]D7p to sum up 
the terms of Eq. ()D6P that are independent of 99;^^, one eventually arrives at Eqs. 

dM]). 



Off-Shell T-matrix of dielectric and magnetic sphere 

Finally, from S{x) and Eq. ()D7p one deduces the expansion 
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which is useful to investigate the limit k ^ kg alluded to in Sec. 14.31 
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